Entanglement renormalization, scale invariance, and quantum criticality 
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The use of entanglement renormalization in the presence of scale invariance is investigated. We 
explain how to compute an accurate approximation of the critical ground state of a lattice model, 
and how to evaluate local observables, correlators and critical exponents. Our results unveil a precise 
connection between the multi-scale entanglement renormalization ansatz (MERA) and conformal 
field theory (CFT) . Given a critical Hamiltonian on the lattice, this connection can be exploited to 
extract most of the conformal data of the CFT that describes the model in the continuum limit. 
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The study of quantum critical phenomena through 
real-space renormalization group (RG) techniques 0,121 
has traditionally been obstructed by the accumula- 
tion, over successive RG transformations, of short-range 
entanglement across block boundaries. Entanglement 
renormalization ^ was recently proposed as a technique 
to address this problem. By removing short-range entan- 
glement at each iteration of the RG transformation, not 
only can arbitrarily large lattice systems be considered, 
but the scale invariance characteristic of critical phenom- 
ena is also seen to be restored 0, |3| • 

In this paper we explain how to use the multi-scale 
entanglement renormalization ansafa (MERA) Q to in- 
vestigate scale invariant systems 0, 0, IH H, 0] • It bas been 
showed that the scale invariant MERA can represent the 
infra-red limit of topologically ordered phases [6]. Here 
we focus instead on its use at quantum criticality. We 
present the following results: (i) given a critical Hamil- 
tonian, an adaptation of the algorithm of Ref. 0] to com- 
pute a scale invariant MERA for its ground state; then, 
starting from a scale invariant MERA, (ii) a procedure 
to identify the scaling operators/dimensions of the the- 
ory and (iii) a closed expression for two-point and three- 
point correlators; (iv) a connection between the MERA 
and conformal field theory, which can be used to readily 
identify the continuum limit of a critical lattice model; 
finally (v) benchmark calculations for the Ising and Potts 
models. 

We note that result (ii) was already discussed by Gio- 
vannetti, Montangero and Fazio in Ref. Q using the bi- 
nary MERA of Ref. [5] . Our derivations are conducted 
instead with the ternary MERA of Ref ^ (see Fig. [T]), 
in terms of which results (iii)-(iv) acquire a simple form. 

We start by considering a finite ID lattice C made of 
N sites, each one described by a vector space V of di- 
mension X- The (ternary) MERA is a tensor network 
that serves as an ansatz for pure states \'^) € V®^ of the 
lattice, see Fig. [TJ Its tensors, known as disentanglers 
and isometrics, are organized in T « log3 N layers, each 
one implementing a RG transformation. Such transfor- 
mations produce a sequence of lattices. 
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where lattice Cr+i is a coarse-graining of lattice Cr, and 
the top lattice Ct is sufficiently small to allow exact nu- 
merical computations. Let a denote a local observable 
supported on two contiguous sites of L, and let pT be 
the density matrix that describes the state of the system 
on two contiguous sites of Ct- Then the ascending and 
descending superoperators At and 2?t- [8], 



AriOr-l), 



Pt-1 = Z'r(Pr), 



(2) 



generate a sequence of operators and density matrices 

A 



Ai A2 
Oq ^ Oi -> 

Po ^ Pi ^ 



ot, 



Go = o, (3) 
pa = P, (4) 



where Or and pr are supported on two contiguous sites 
of the lattice Ct- Eq. ^ allows us to monitor how the 
local operator o transforms under successive RG trans- 
formations, whereas its expected value (o) — tr(po) can 
be evaluated by computing p in Eq. (|4]). 

RG fixed point. — The scale invariant MERA cor- 
responds to the limit of infinitely many layers, T — > 00, 
and to choosing the disentanglers and isometries in all 
layers to be copies of a unique pair u and w [1, Q . In this 
case we refer to the ascending superoperator At, which 
no longer depends on r, as the scaling superoperator S 
(see Fig. [T]), and to its dual Vt as S* - Notice that 5 is a 
fix ed-p oint RG map. Then, as customary in RG analysis 
@, [3 1 the scaling operators (pa and scaling dimensions 
Aa of the theory, 

S{(f>a) ^ Xa(f>a, A^^-loggAa, (5) 

are obtained by diagonalizing this map. 



5(.) = ^A„ 



»), tr{(f>a(l)f3) = 5ap, (6) 



where (j)a are the eigenvectors of the dual S* , S*{<j)a) — 
Xa<j>a- Eq- [Hwas first discussed__in Ref. 3 by Giovan- 
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Ci 



Ct, 



Co^C, (1) 



netti, Montangero and Fazio 
vious observation (see Eq. 5 of Ref. 



It formalizes a pre- 
0]) that the scale 
invariant MERA displays polynomial correlations. By 
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(iii) 




FIG. 1: (Color online) (i) Two lowest rows of disentanglers u 
and isometries ui of the ternary MERA. They map the orig- 
inal infinite lattice Co = jC into increasingly coarse-grained 
lattices jCi and £2- Notice that three sites of £t-i become 
one site of £t, hence the use of log3 throughout the paper. 
{ii)-{iv) Under the coarse- graining transformation defined by 
the MERA, two-site operators supported on three different 
pairs of sites of £t-i become supported on the same pair of 
sites of £t- (v) Accordingly, the scaling superoperator S is 
the average of three contributions, each of which (and thus 
also their average) is unital and contractive thanks to the 
isometric character of u and w [51. 



construction, S is unital, iS(I) 



so that the identity 



operator I m 



is a scaling operator with eigenvalue 



Ai = 1; and contractive, meaning |Aq| < 1 |12l |. Here we 
will assume, as it is the case in the examples below, that 
only the identity operator I has eigenvalue A = 1. Then 
the operator p = I is a density matrix that corresponds 
to the unique fixed point of S* , S* (p) = p, and since 



lim ( S* o- - -oS\ ){pt) 

T times 



p 



for any starting pr, it follows that p is the state of any 
pair of contiguous sites of C. [Consistent with scale in- 
variance, p is also the state of any pair of contiguous sites 
of Cr for any finite r] . The computation of the expected 
value of the local observable o is then straightforward. 



(o) = tr(/5o), 



(8) 



which for the scaling operators reduces to {4>a) — 5ai- 

Correlators. — Let us now diagonalize the one-site 
scaling superoperator S^^^ of Fig. [2l 



(9) 



where the scaling dimensions IsSa'^ = — logg Aa ^ coincide 
with Aq |T3|. The correlator for two scaling operators 




o'"(.v) o'"(v-) 



FIG. 2: (Color online) (i) One-site operators on special sites 
are coarse-grained into one-site operators, (ii) Scaling super- 
operator for one-site operators, (iii) In computing correlators 
on specific sites x and y (or x, y and z), one-site operators 
are coarse-grained individually according to 5'^' until they 
become nearest neighbors (which in this case occurs at lattice 
^2, 9 = 2). 



and (j)"}'^ placed on contiguous sites reads 



C^p ^ (</.W(l)0(^'(O)) = tr((4i) ® (10) 



(1) 



Suppose now that 0q and (1)^^' are placed in two spe- 
cial sites X, y as in Fig. [21 where r^y = x — y is such that 
\fxy \ = 3' for (7 = 1, 2, • • • . Then after q = logg \rxy \ itera- 
tions of the RG transformation, and (jJ^'^ become first 
neighbors again. Notice that each iteration contributes a 
Using the identity a'°s'' = we find 



factor Aq''A^^\ 



l-A(,i)-A<,^' 



and obtain a closed expression for two-point correlators 

Cap 



(7) For three-point correlators we define the constants 



"a/3 



Pi 



.(1) 



A(i) 



(11) 



(12) 



tr((0«0(/.«®0«)p(3)) (13) 



where the trace corresponds to the correlator on three 
consecutive sites and 

p(3) is obtained from p. For \rxy \ — 
''xz\l'^ = 3'', analogous manipulations lead to 



o ° I 

yz\ I ^zx 



-(14) 



CFT. — The continuous limit of a quantum criticial 
lattice system (scale invariant case) corresponds to a con- 
formal field theory (CFT) A CFT contains an in- 
finite set of quasi-primary fields (j)'^'^ , with scaling dimen- 
sions A^'^'^. The correlators involving two or three quasi- 
primary fields have expressions analogous to Eqs. Illl and 
[141 the (symmetric) coefficients C^p^ for three-point 



3 



correlators coincide with those in the so-called operator 
product expansion (OPE). Moreover, quasi-primary fields 
are organized in conformal towers corresponding to irre- 
ducible representations of the Virasoro algebra. Each 
tower contains one primary field at the top, with 
conformal dimensions {t, i) [such that its scaling dimen- 
sion is A^* = t + t], and its infinitely many descendants, 
which are quasi-primary fields with scaling dimension 
A = AP -t- n for some integer n > 1. 

A CFT is completely specified by its symmetries once 
the following conformal data has been provided: (i) the 
central charge c, (ii) a complete list of primary fields with 
their conformal dimensions and (iii) the OPE for these 
primary fields. For instance, the Ising CFT in 1+1 di- 
mensions has central charge c — 1/2, three primary fields 
identity I, spin a and energy e with conformal dimensions 
(0,0), (j^, j^) and (i, i), and OPE coefhcients 



i^CFT c ^CFT 



2 ' ^crcrcr ' 



CFT , 

eee 



CFT _ 

eecT 



0. (15) 



The present analysis readily suggests a correspondence 
between the scaling operators 0q of the scale invari- 
ant MERA, defined on a lattice, and the quasi-primary 
fields (j)'^^'^ of a CFT, defined in the continuum. To- 
gether with the algorithm described below, this corre- 
spondence grants us numerical access, given a critical 
Hamiltonian H on the lattice, to most of the conformal 
data of the underlying CFT, namely to scaling dimen- 
sions and OPE coefficients. The central charge c can 
also be obtained e.g. [l3| from the von Neumann en- 
tropy S{p) = — tr(plog2p), which for a block of L sites 
scales, up to some additive constant, as S* = | log2 L p^ . 
We then have S{p) - S{p'^^'>) = f (logj 2 - logj 1) = f , or 
simply 



c = 3 



(5(p)-5(/5«) 



(16) 



Algorithm. — Given a critical Hamiltonian H for an 
infinite lattice, we obtain a scale invariant MERA for 
its ground state by adapting the general strategy 
discussed in Ref. [1]. Recall that tensors (disentanglers 
u and isometrics w) are optimized so as to minimize the 
energy E = {'i>\H\'i>). After linearization this reads 



E — tr{uTu) + ki = tr(wT„) + k2, 



(17) 



where T„ and Ty, are known as environments and fci , k2 
are two irrelevant constants. In the translation invariant 
case Q the environment for, say, an isometry w at layer t 
of the MERA, — f{ur, w^-, p^, /i^_i), is a function of 
the disentangler Ut and isometry Wr of that layer, a two- 
site density matrix p-r and a two-site Hamiltonian term 
/it-1- In the present case, we replace the above with 
the unique pair {u,w), the fixed-point density matrix p, 
and an average Hamiltonian h = h.,- /i'^ , where the 
weights 1 /y account for the relative number of tensors in 
different layers of the MERA. Then, starting from some 
initial pair (u, w) and the critical Hamiltonian H made of 
two-body terms h, the following steps are repeated until 
convergence: 
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FIG. 3: (Color online) Scaling dimensions obtained from 
the spectrum of the scaling superoperator <S. Circles indicate 
primary fields. Left: For the Ising model we can identify the 
scaling dimensions of the three primary fields, the so-called 
identity I, spin a and energy e, together with several of their 
descendants. Right: The spectrum of S for the 3-level Potts 
model shows some of its primary fields, including its primary 
fields with multiplicity two, namely the spins ai and (T2 and 
the pair Zi and Z2 [10l |. 



Al. Given the latest {u,w), compute {p,h). 

A2. Given {u,w, p,h), update the pair {u,w). 

In step Al, the scaling superoperator S is built as in- 
dicated in Fig. [T] We compute the fixed-point density 
matrix p by sparse diagonalization of S, and the average 
Hamiltonian h by using hr — S{hr-i), Hq = h [T6|. Step 
A2 is decomposed into a sequence of alternating opti- 
mizations for u and w as in the generic algorithm of Ref. 
Q , where each tensor is updated by computing a singular 
value decomposition of its environment. 

Examples. — We illustrate the above ideas and the 
performance of the algorithm by considering the Ising 
and 3-level Potts quantum critical models in ID, 



'■1 ' crM (!!'■+ 1 



Z ' ^ X ^ X 



(18) 



Potts 



where cr^ and ax are Pauli matrices, and 
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Mx.2 = {Mx.i)'^- Notice that sites have a vector space of 
dimension c? = 2 or d = 3. In order to use a scale invari- 
ant MERA with x > d, we allow the disentanglers and 
isometrics of the first few (typically one to five) layers to 
be different from u and w. We iterate steps A1-A2 about 
1000 times. With a cost per iteration that scales as 
and using a 3 GHz dual core desktop with 8 Gb of RAM, 
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Ising 


^CFT 


A (MERA X = 22) 


rel. error 


o 


1/8 = 0.125 


0.124997 


0.002% 


e 


1 


1.0001 


0.01% 


Potts 


^CFT 


A (MERA X = 22) 


rel. error 


o"i 


2/15 = 0.13 


0.1339 


0.4% 




2/15 = 0.13 


0.1339 


0.4% 


e 


4/5 = 0.8 


0.8204 


2.5% 


^1 


4/3 = 1.3 


1.3346 


0.1% 


^2 


4/3 = 1.3 


1.3351 


0.1% 



TABLE L Comparison of scaling dimensions of primary 
fields of the Ising and Potts models calculated using MERA 
(A(MERA X = 22)) with exact results known from GET 
(A<=""^). 



simulations for x = 4, 8, 16, 22 take of the order of min- 
utes, hours, days and weeks respectively. The following 
results correspond to % = 22. 

From Eq. [TB] we obtain an estimate for the central 
charge, namely cising = .5007 and Cpotts = -806, to be 
compared with the exact results 0.5 and 0.8. Fig. [3]shows 
the smallest scaling dimensions of the scaling super- 
operator S [l^l- We obtain remarkable agreement with 
those expected from CFT, as shown in Table HI Recall 
that all the critical exponents of the model can be ob- 
tained from the scaling dimensions of primary fields. For 
instance, for the Ising model the exponents v and 77 are 
V — 2Ao- and 77 — 53^1 whereas the scaling laws express 
the critical exponents a, (3,^,5 in terms of v and 77 p^ . 
Further, the OPE coefficients for primary fields of, say, 
the critical Ising model are computed as follows. The ma- 
trix Caf3 in Eq. [10] is diagonal for the scaling operators 



corresponding to I, a and e, which we normalize so that 
Ca/3 = ^ap- With this normalization, we then compute 
the coefficients Cap-y using Eq. [131 We reproduce all the 
values of Eq. [T5|with errors bounded by 3 x 10~^. 

Discussion. — In this paper we have explained how 
to compute the ground state of a critical Hamiltonian us- 
ing the scale invariant MERA and how to extract from it 
the properties that characterize the system at a quantum 
critical point. Our results, which build upon those of Ref. 
[E 0, SH, 0, [1| , also unveil a concise connection between 
the scale invariant MERA and CFT. This correspondence 
adds significantly to the conceptual foundations of en- 
tanglement renormalization. The scale invariant MERA 
can be regarded as approximately realizing an infinite di- 
mensional representation of the Virasoro algebra joj, ■ 
The finite value of x effectively implies that only a finite 
number of the quasi-primary fields of the theory can be 
included in the description. Fields with small scaling di- 
mension, such as primary fields, are retained foremost. 
As a result, given a Hamiltonian on an infinite lattice, 
we can numerically evaluate the scaling dimensions and 
OPE of the primary fields of the CFT that describes the 
continuum limit of the model. This approach differs in 
a fundamental way from, and offer an alternative to, the 
long-established techniques of Refs. based instead 

on finite size scaling. 

We conclude by noting that most of our considerations 
rely on scale invariance alone and can be applied to study 
also critical ground states in 2D systems [18| . 
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